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Abstract 

Many theories beyond the standard model (BSM) contain new CP-odd and CP-even neutral 
scalars 4> = {^5 and new vector-like fermions The couplings of the CP-odd scalar A to two 

standard model (SM) gauge bosons cannot occur from renormalizable operators in a CP-conserving 
sector, but can be induced at the quantum loop level. We compute these effective couplings at the 
1-loop level induced by the SM fermions and vector-like fermions, present analytical expressions for 
them, and plot them numerically. Using the 8 TeV Large Hadron Collider (LHC) 77 , and 

ti channel data, we derive constraints on the effective couplings of the cj) to standard model gauge 
bosons and fermions. We present the gluon-fusion channel cross-sections of the 0 at the 8 and 14 TeV 
LHC, and its branching-ratios into SM fermion and gauge-boson pairs. We present our results first 
mo del-independently, and then also for some simple models containing cj) and 'ipvL in the singlet and 
doublet representations of SU{2). In the doublet case, we focus on the two-Higgs-doublet (2HDM) 
Type-H and Type-X models in the alignment limit. 


1 Introduction 

A long series of experiments culminating in the Large Hadron Collider (LHC) discovery of the Higgs boson 
at a mass of about 125 GeV has firmly established the standard model (SM) as the correct description 
of Nature up to an energy scale of a few hundred GeV. With this discovery, the theoretical puzzle as 
to why the Higgs boson remains this light when quantum effects should correct it to the highest scales 
present in the theory (such as the Planck scale) comes to the fore. This problem of the stability of the 
electroweak (EW) scale is the well known hierarchy problem of the SM. This could be a clue that some 
new physics beyond the standard model (BSM) is present near the EW scale which renders it stable 
against quantum corrections, making it natural. Many theoretical proposals have been made for this 
new physics (for reviews see Ref. [1]), and they usually contain new particles at the TeV energy scale. 
We are poised at a very interesting time when the LHG is probing this energy scale and can tell us if 
one of these proposals is realized in Nature. 

Among the possibilities of BSM physics that makes the EW scale natural are models in which the 
Higgs-doublet of the SM is a pseudo-Nambu-Goldstone boson (pNGB). Goncrete realizations of this 
idea, for example, are in models of little-Higgs, composite-Higgs and extra dimensions (for reviews see 
Refs. [2, 3, 4] respectively). In such models, in addition to the GP-even Higgs boson, there could be 
new GP-odd scalar (A) and GP-even scalar (iL), which we denote collectively as 0 = {A, iL}, that are 
also pNGBs due to which their mass is much lower than the cut-off scale. Also, new heavy vector-like 
fermions (VLE, denoted as ^vl) are usually required, that along with the SM fermions (SME), complete 
some representation of a bigger group containing SU (2) (8)17(1). The new vector-like fermions can include 
vector-like quarks (VLQ) and vector-like leptons (VLL) and may be present in addition to the usual SM 
quarks (SMQ) and leptons (SML). By vector-like fermions we mean that fermions in a representation of 
the SM gauge-group and in its conjugate representation both appear in the theory (for more details see 
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for example Ref. [5]). Some supersymmetric models also include vector-like matter, and thus have (j) and 
2pvL both present, along with many superpartners. 

The phenomenology of a CP-odd scalar at the LHC can be quite distinct as compared to a CP-even 
scalar (such as the SM Higgs boson), and one focus of this work is to elucidate this aspect. If CP- 
invariance is not spontaneously broken by an A vacuum expectation value (VEV), i.e. if (H) = 0, as 
we assume here, AW~^W~, AZZ (collectively called AVV couplings), and also H 77 and AZ^ couplings 
cannot arise from renormalizable operators. Also, the last two do not arise from renormalizable operators 
because of unbroken electromagnetic (EM) gauge invariance, the same reason why ^77 and hjZ are zero 
at the renormalizable level. These can then only result from higher-dimensional operators generated at 
loop-level. In contrast, for the CP-even SM Higgs boson (denoted as h), the hW^W~ and hZZ couplings 
are generated at tree-level from dimension-four operators after electroweak symmetry breaking (EWSB), 
i.e. with (h) = n/\/2. Therefore, generically speaking, the AW~^W~ and AZZ effective couplings, 
generated at loop-level, are much smaller in magnitude compared to the tree-level hW^W~ and hZZ 
couplings; the A 77 and ^77 effective couplings are both loop suppressed and small, and similarly the 
AjZ and hjZ are also both loop suppressed. Thus, similar to the h, the gg ^ A “gluon-fusion” 
channel is important at the LHC, while compared to the h, the vector-boson fusion channel of A is much 
suppressed. The alternate possibility of (A) 7 ^ 0 is not discussed here but is considered for instance in 
Refs. [ 6 ]. 

Turning next to 0 = {A, H} couplings to fermions, we include (j) couplings to new vector-like fermions 
at the tree-level. Eurthermore, if 0 is part of a doublet, it couples also to SM fermions at the tree-level 
(similar to h). We consider the case when (j) couples significantly only to third generation SM fermions, 
a situation common in many BSM extensions. Thus, the relevant couplings to SM fermions are 066, 
and If the 066 coupling is sizable, 66 (j)^ bg ^ 60 and gg bbcj) can be important 

production channels of the 0. However, we do not include these production channels in this work, but 
restrict ourselves only to the gluon-fusion channel. 

Mostly, we restrict ourselves to the situation when < 2Myl so that 0 cannot decay to a pair of 
VLEs. If the 'ipvL is light enough they can also be studied directly at the LHC as discussed for instance 
in Ref. [5] and references therein. However, if they are too heavy to be directly produced at the LHC, 
but the 0 (or h as studied in Ref. [7]) can be directly produced and it’s couplings measured, the VLE 
contributions to the 0 effective couplings we derive here can be useful in probing the ^pvL indirectly. 

We identify the lighter CP-even state (h) to be the 125 GeV state discovered, and whose properties 
measured, at the LHC. The h couplings measured at the LHC so far largely agree with the SM, at 
least to about a few tens of percent, and the magnitude of the hVV coupling (with V = is 

constrained to be close to the SM coupling at the few tens of percent level. This will be realized in the 
so called “decoupling limit” [ 8 ] , or more generally in the “alignment limit” [9] . In order to capture many 
different BSM models, we perform a model-independent effective theory analysis of the 0 coupled to SM 
fields. We present the constraints from the recent 8 TeV LHC run using the 77 , and ti channels, 

and present the signal cross section (c.s., a) at the LHC as a function of the effective-couplings of the 
0 (denoted by hz) and the branching ratio {BR) into these modes. We do not focus much on the ZZ 
and W^W~ decay channels of the 0 as the branching-ratio into these modes are much smaller than the 
other modes due to AVV coupling being generated only at the loop-level, and the HVV coupling being 
zero in the alignment limit. We also present many simple models containing A and 'ipvL in SU( 2 )-singlet 
and doublet representations. Eor A in a doublet, we restrict ourselves to the two-Higgs-doublet model 
(2HDM) Type-H and Type-X. We present the 1-loop analytical expressions for the {A^f^f, A 77 , A 7 Z} 
effective couplings induced by SMEs and VLEs in each of these models; as a function of the model 
parameters, we plot numerically these effective couplings and the BR into the 77 , jZ and fermion final 
states. These are some of the main results of this work. 

In previous studies, one of us has considered the implications of models with VLQs and VLLs coupled 
to the lighter CP-even Higgs boson h in Ref. [7], and the direct LHC signatures of VLQs in Refs. [5]; this 
work complements them by considering aspects of heavier neutral CP-odd and CP-even scalars A, H. 
In Ref. [10] we study many aspects dealt with in this paper but in a specific little-Higgs model, the 
SU( 6 )/Sp( 6 ) model by Low, Skiba and Smith [ 11 ]. We also list there many little-Higgs models that 
contain a 2HDM structure. The results of this paper are useful in deriving constraints and prospects of 
such models. 

Erom the vast literature, we give a sampling below of studies that deal with extra BSM neutral 
scalars, have overlap with our work and that take into account the recent LHC 8 TeV constraints. We 
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also mention how our work complements them. There exist several studies which present a{pp A) (see 
for example Refs. [ 12 , 13]) in the context of 2 HDM, minimal supersymmetric standard model (MSSM) 
and next-to-MSSM (NMSSM). We highlight the effects of VLFs on a{gg A) in various SM extensions 
including 2 HDM-II and 2 HDM-X. Refs. [14, 15] consider the possibility that the observed 125 GeV state 
at the LHC is a CP-odd scalar, and the former shows that this possibility is disfavored by the LHC data. 
Refs. [16, 17] analyze 2 HDM Types I and H taking into account the 125 GeV LHC data, all pre-LHC 
constraints and results of the heavy-Higgs searches in various channels. Ref. [18] performs a global fit 
of general 2 HDMs using ATLAS, CMS and Tevatron results. Ref. [19, 20, 21 , 22 , 23] shows the allowed 
parameter space of 2HDM-H, applying theoretical (perturbativity, unitarity and vacuum stability) and 
experimental (LEP, Tevatron and LHC 125 GeV Higgs data, precision observables and H-physics and 
electric dipole moment measurements) constraints. Ref. [24] also includes the heavy Higgs exclusion 
limits to constrain the 2 HDM. LHC 8 TeV constraints on the 2 HDM parameter-space are also discussed 
in Refs. [25, 26, 27, 28, 29]. The heavy neutral scalars of the 2 HDM, namely A and iL, are studied in 
Ref. [30], where the LHC 8 TeV exclusion and 14 TeV reach from the processes gg ^ H ^ AZ and 
gg ^ A ^ HZ are presented. Ref. [31] constructed an 5'0(5) symmetric 2HDM which naturally realizes 
the ’’alignment limit” and puts constraints on it’s parameter space from the 8 TeV LHC data. Ref. [32] 
puts limits on the the triple Higgs couplings and presents a set of benchmark points for probing SM-Higgs 
pair production and the search of heavy Higgs bosons through non-standard decay channels (i.e decays 
of A, H that involves at least one Higgs boson in the final state). Ref. [33] calculates the loop factors for 
the AVV couplings in the MSSM and the 2 HDM with a heavy chiral fourth generation. Ref. [34] studies 
A WW^ ZZ decays and compares this with the corresponding CP-even scalar decays in 2 HDM-H, and 
also with a chiral fourth generation or additional heavy vector-like quarks (VLQ) added. In addition to 
these, here we also include the effects of VLFs on A ^ 77 , Z 7 decays. An effective Lagrangian analysis 
of new heavy scalar particles is presented in Ref. [35] . Various VLF models and related phenomenological 
issues are also studied in Refs. [36]. Many of these studies are done with specific models in mind while we 
present the LHC limits and signal c.s. in a model-independent manner, and using these, derive results 
for the models we introduce, and also for some of the models above. 

The paper is organized as follows: In Sec. 2 we present a model-independent analysis of the CP-odd 
and CP-even neutral scalars 0, present constraints on its effective couplings from the 8 TeV LHC run, 
the LHC gluon-fusion c.s., and BR into SM fermion and gauge boson decay modes. In Sec. 3 we present 
many simple models containing (j) and ipvL as SU(2) singlets or doublets. For each of these models, 
we work out the 1-loop effective couplings of the (j) and present its BR into two body decay modes. 
One can read-out the current constraints and gluon-fusion c.s of the (j) at the LHC for each of these 
models in conjunction with the results in Sec. 2 . The models considered include 0 as an SU{2) singlet, 
or contained in the 2HDM, with correspondingly the 'ipvL also in singlet or doublet representations. We 
offer our conclusions in Sec. 4. For the various models we discuss, we compile expressions for the mass 
eigenvalues and mixing angles in App. A, and the 1 -loop effective couplings in App. B. 

2 Mo del-independent Analysis 

In this section, we define an effective Lagrangian with couplings of the neutral scalars, CP-odd A and CP- 
even h, H to SM gauge bosons and fermions. We denote the neutral scalars collectively as (j). In models 
that contain two CP-even scalars, we identify the lighter one (h) as the 125 GeV scalar observed at the 
LHC. For the heavier states {A^H), we show the constraints from the 8 TeV LHC, signal c.s. ax BR 
into various SM two body final states at the 8 and 14 TeV LHC, as a function of the effective couplings 
and rriff). For any given new physics model, one can obtain this effective Lagrangian by integrating out 
heavier fields, following which the results of this section can then be used to obtain the LHC limits and 
gluon-fusion cross-section in that model. 

CP invariance requires the CP-odd scalar A coupling to SM gauge bosons to be only via higher 
dimensional operators. The CP-even scalars can couple to the massive gauge bosons at tree level. 
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( 1 ) 


Showing only the new physics terms, the effective Lagrangian for any neutral scalar 0 is 

+ y^WwWW^ + y^PZZCPZ^^Z^ - ^^^K^-y-y4>Y^uarF^^F^’' - K^-yZ4>Y^uarF^^Z^’' 

where X = 75 ,F^zycrr = ^^xuar for the CP-odd scalar, while X = I (identity matrix), Y^^ctt = ducrdvr 
for the CP-even scalar. Here i^cf^ij s contain other fermion and gauge boson loop contributions. Tree 
level scalar gauge boson couplings y(j)ZZ^y(j)WW are zero for the A. We have defined the dimensionless 
effective couplings k by pulling out a new-physics mass-scale M in the effective (pVV terms. For the 
numerical results we show, we set M = 1 TeV from now on and show only /^, and for other values of M, 
the K can easily be rescaled. Although we have defined the effective couplings k. by extracting a heavy 
new-physics mass scale M, SM fermion contributions are to be included when present. Eq. ( 1 ) is an 
effective Lagrangian at a scale just above m^. Heavy BSM fermion and the SM fermion contributions 
are to be included in k before comparing with the plots we show in this section. For various simple SM 
extensions detailed in Sec. 3 we compute the /^’s and present them in App. A. If SM fermions contribute 
and can go onshell, the k are complex. In this case, the Kcpvv that appear in our plots in this section 
should be read as \i^(f)Vv\- We assume to be real in this work. 

The CP-odd scalar can decay to SM gauge bosons or fermions. In terms of the hz and y^s defined 
above, the decay rates to different final states are 

r(^ ^ Z7) = (iSSf)' ^ ^ 

where n = 3 and n = 1 for CP-even and CP-odd scalars respectively, Vf = rz = with 

Nc = 3 for quarks and 1 for leptons. Here we have defined r(0 ^ gg) to have an extra factor of 8 
compared to r(0 ^ 77 ) anticipating a color factor. It turns out however that for a quark in the loop, 
the color factor in the r{(j) ^ gg) is actually 2. This will get compensated for in hz^j^gg (see for example 
Eq. (B.34)). Using these expressions, one can work out the BR of the 0 into these final states in any 
new physics model. 

We turn next to discussing limits from the 8 TeV LHC and the gluon-fusion cross-section at 14 TeV. 
To obtain the limits on the effective couplings hz and we use upper-limits {UL) from recent LHC analysis 
on (7{pp (j))x BR {(j) XX), and the currently relevant constraints are XX = { 77 ,T+r“,tt}. We 

take the limits on the 77 channel from the CMS analysis Ref. [37] which has an upper limit up to of 
850 GeV, on the r+r“ channel from the ATLAS analysis Ref. [38] up to of 1000 GeV, and from the 
ATLAS analysis Ref. [39] for the tt channel. Using these we constrain the effective couplings of Eq. ( 1 ). 

At the LHC, the 0 can be produced hy gg ^ (j) (called gluon-fusion channel), which starts at the 
1-loop level when cj) couples to colored fermions. In addition to the above production channel, if 0 couples 
to 5-quarks, there are additional production channels, namely, bh ^ (j) (called 55-fusion), hg bcj) and 
gg bbcj) (called 5-quark associated production) channels; how these compare with the gluon-fusion 
channel depends on how large the 550 coupling is in a given model. For instance, for = 0.5, we find 
that the production rate via 55-fusion and 5-quark associated production channels becomes comparable 
to the gluon-fusion channel with tz^^gg ^ 20. We include only the gluon-fusion channel in this study, 
but in models with a large 550 coupling, the 55 fusion and 5-quark associated production channels may 
have to be included, which we do not do here. For a study involving the 5-quark associated production 
channels of the h including gg bbh, see Ref. [40]. One can separately study the 5-quark associated 
production channels by tagging on the final state 5-jet as discussed in Ref. [38]. Ref. [41] has recently 
studied 55 fusion and 5-quark associated production channels for a light CP-odd scalar. Although there 
are some LHC limits using 5-tagged events to which the 55 decay mode and the 5-quark associated 
production channels contribute, we do not include them in our analysis here. So far these results have 
been presented for < 350 GeV (see Refs. [42, 43, 44]). 

Rather than compute the A, H production rate at the LHC ourselves, we relate it to the SM Higgs 
production rate at the same mass, and make use of the vast literature on h production rate. Since 
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Figure 1 : cr{gg (j)) (in pb) at the 8 TeV LHC (left) and 14 TeV LHC (right) for = 200 GeV (red), 
500 GeV (blue), 800 GeV (green) and 1000 GeV (yellow). 
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Figure 2 : 8 TeV LHG constraints from the 77 channel (left), channel (middle) and tt channel 

(right), for = 200 GeV (red), 500 GeV (green), 800 GeV (blue) and 1000 GeV (yellow). The regions 
to the top and right of the curves are excluded at the 95 % GL level. 


cr{gg -^ (/)) oc T{(/) ^ gg)^ we can write the a * BR for (j) production followed by decay into the final-state 
W as 

a{gg (j)) = ^ x <^{99 ^ h) . (3) 

From Eq. (2) we compute F((/) ^ gg) for a given and using Eq. (3) we apply the upper-limit (UL) 
from the 8 TeV LHC quoted above for various BR{(j) XX). Eor our numerical work, we calculate 
T{h gg) in the SM and take cr{gg h) from Ref. [45]. We assume here that the dependence on the 
PDE, and the acceptance at the LHC for H, iL and h are not very different, which should be reasonable 
assumptions. Eor the decay A ^ XX^ the final-states W we consider are 77 , and ti as these are 

currently the significant ones. We compute the BR{A XX) using Eq. (2). If A^H are fairly close in 
mass, i.e. closer than the experimental resolution to separate them (say 30% of m^), and no kinematic 
variables can separate them, we should include all of them into the a * BR above. 

In Eig. 1 we show cF{gg (j)) at the 8 TeV LHC (left plot) and 14 TeV LHC (right plot) as a function 
of i^(f)gg- o-{gg (j)) is obtained using Eq. (3) and the (j{gg h) from Ref. [45] as mentioned earlier. In a 
given new physics model, one can compute Hi^gg and then use these plots to obtain the a{gg (j)). Using 
the a{gg ^ 0), we obtain constraints from the 8 TeV LHC data as a function of the BR into a particular 
mode. We show this in Eig. 2 obtained from the 77 , and tt channels. The regions to the top 

and right of the curves are excluded at the 95 % CL level. In the 77 channel, the bound is strongest for 
rricf) = 200 GeV since the experimental exclusion is tightest at that mass. We see that there is no constraint 
from this channel for BR{(j) 77 ) < 10“^ for the range of n^gg shown. Erom the channel, we find 

the strongest limit for of about 500 GeV since the experimental exclusion is tightest at that mass. 
We show in Eig. 3 the total cr{gg ^ 0) * BR{(j) XX) contours (in pb) for XX = { 77 , r+r“, tt} at the 
14 TeV LHC, making use of the fact that the total cr{gg (j) ^ XX) oc * BR{(f) XX), omitting 
kinematic factors independent of couplings. Thus, each mode XX can be considered and presented 
independently of the others as we do here. The 95 % CL LHC exclusion discussed above is also shown 
labeled as ’8 TeV’. If the (pbb coupling is large, i.e. bigger than about 0.5, inclusion of the 6 -fusion and 
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6-associated production channels (along with the (pgg channel that we have included here) could result 
in a stronger exclusion than we obtain here. 

As already mentioned, the mo del-independent results presented in this section can be used to ob¬ 
tain the LHC constraints and gluon-fusion c.s. in any particular model by computing first the effective 
couplings in that model. We compute the effective couplings in many simple models next. 

3 Models 

In this section we consider some specific models for the neutral CP-odd and CP-even scalars A, H and 
study their LHC production and decays into 2-body final states. We compute the decay rates assuming 
a sharp turn-on at threshold of the 2-body final state. The goal is to capture in simple models many 
of the features present in realistic BSM models as far as the LHC phenomenology of A, H is concerned. 
As before, we collectively denote A^H as (j). We mostly focus on the situation when < 2Myl and 
do not focus on the phenomenology due to the (j) decaying to a pair of on-shell VLF. We first consider 
the models where 0 is an SU (2) singlet and couples to SU (2) singlet VLF (singlet A with a vector-like 
up-type-singlet or SVU model, or with a down-type-singlet or SVD model) and SU{2) doublet VLF 
(singlet A with minimal vector-like quark doublet or SVQ model). We next consider effective models 
with (j) in an SU{2) doublet, with the two SU(2) doublet scalars and ^2 both having hypercharge 
+1/2. The 2HDMs we consider are either Type-H like or Type-X like. We notate the Type-H like 
models, for example, as MVQD for minimal vector-like extension with VLQ doublet Q and down-type 
VLQ singlet D, and MVQU for a similar model with an up-type VLQ singlet U instead, and a similar 
model with the 2HDM Type-X structure instead as MVQDX. We include subscripts depending on 
which Higgs doublets the fermions couple to, i.e. MVQDij will mean that the model has one VL-quark 
doublet 7 /, and one down-type VL-quark singlet y, with the couplings 'ipLXn^i turned on. 

Among our example models are some that mimic BSM models that have (p Yukawa couplings with an 
SMQ and a VLQ, for example, the 3rd generation SMQ with an up-type singlet-VLQ to give the MVU 
model. 

Many of the effects we present are similar for the CP-odd and CP-even scalars A, H. One important 
difference between the A and H is that at tree-level, the AVV (with V = {W, Z}) couplings are zero and 
are only generated by SM and BSM fermions at the loop level, while the HVV couplings could be nonzero 
at tree-level. However, in the alignment limit we consider (discussed later), the HVV couplings are zero. 
Although the V{H ^ 77 ) gets a contribution from the charged scalar {H^) loop, while r(A ^ 77 ) does 
not, this is very small compared to the fermionic contributions [46]. Thus in the alignment limit the A 
and H have very similar phenomenology. Therefore we will mostly present the phenomenology of the 
CP-odd scalar A, and where relevant, we will also contrast it with the situation for the H. Since the 
tree-level HVV coupling is zero in the alignment limit, in the appendix we only give the expressions for 
the fermion contributions to the For the SM Higgs we must include the W loop contribution to 

and RhZ'y which we do not present here. 

3.1 Model with an SU{2) singlet A with VLQ-VLQ Yukawa couplings 

We start by considering some models with an SU ( 2 ) singlet A coupled to SU ( 2 ) singlet or doublet VLFs. 
For an SU( 2 ) singlet (j) one cannot write Yukawa couplings with chiral SMFs, and thus ggcj) and 77 ^ 
couplings can only be induced by VLFs, if they are present, as we explicitly show here. Thus, LHC 
signals of a BSM singlet 0 becomes possible if colored VLFs are coupled to it. 

SVU model : We study a model, which we call 5'W-model, with an SU{2) singlet CP-odd scalar A, 
coupled to an SU{2) singlet, SU{3) triplet VLQ ( 7 /) with hypercharge Clearly, the electromagnetic 
charge Q = Y^. To the SM Lagrangian we add 

£ + eQAf.'ipYi’ - 9Q—Z^'ipY'>P (4) 

Z Z Cw 

+ 'ijjilp'il; — iyAA^jj'y^^p — — ^^A^H^H. 

6 

^ A model with only a vector-like lepton singlet is uninteresting for A phenomenology since no LHC production channels 
are significant (note that the Abb coupling is also not possible in this case). 
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Figure 3: Contours of the 14 TeV LHC a x BR (in pb) in the 77 channel (upper-row) and the 
channel (middle-row) for = 200 GeV (left), 500 GeV (middle), 800 GeV (right), and in the tt channel 
(bottom-row) for ttia = 500 GeV (left), 800 GeV (middle), 1000 GeV (right). The region to the right of 
the contour labeled ’8 TeV’ is excluded at the 95 % CL level from 8 TeV LHC result. 
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Figure 4: BR (^4 —> 77 ) (black), BR (^4 ^ 7 Z) (blue), BR(74 ^ ZZ) (red), BR(74 ^ 14^14^) (cyan) as a 
function of tra with ha = 0.1 and = 1000 GeV for SVU (left) and SVQ (right) models. 



m^(GeV) 

Figure 5: RAggly\ as a function of tyia for = 800 GeV (red) and 1200 GeV (blue) for SVU model. 


The SM Higgs doublet is written as H here. Here we have not considered possible terms coupling the 
H to a SM fermion and a VLF for = 2/3,—1/3 such as qHipR. We study this 

possibility of off-diagonal couplings between the 3rd generation SMQ and a VLQ in the context of the 
SU( 2 ) doublet 4> in Sec. 3.2.3. 

We restrict ourselves to rriA < 2 Myl, so that A cannot decay to a VLF pair. The possible decay 
modes of A are to gg, 77, Z7 and ZZ through a VLF loop, but no decay to W^W~. A cannot decay 
to a pair of SM fermions since such couplings are forbidden by gauge invariance. The effective AV^V^ 
couplings induced by VLFs are given in App. B. From these we compute the partial widths and the BR 
into the above modes. In Fig. 4 we plot BR(A ^ 77), BR(A ^ Z7) and BR(A ^ ZZ) where we chose 
V/; = 2/3 as an example. BR(A ^ gg) is almost constant at around 0.999. 

In Fig. 5 we plot KAgglv\ as a function of From this, one can read-off the (j{gg A) at the 
8 and 14 TeV LHG from Fig. 1 in Sec. 2 . The peaks in Fig. 5 are due to the VLFs going onshell, although 
as mentioned earlier, we do not explore its consequences in this work. In this model, the gluon-fusion 
c.s. of A is induced only through loops of the heavy VLFs due to which the 8 TeV LHG exclusion limits 
on cr X BR into the ZZ channel (see Ref. [47]) or the 77 channel (see Ref. [37]) are rather weak, unless 
Pa becomes so large that perturbativity is lost. 

If tua < 'rnhl‘^ (where h is the 125 GeV Higgs), then h AA becomes kinematically allowed and 
becomes a means of producing A in addition to the gluon-fusion channel discussed above. In Fig. 6 we 
plot BR(h ^ AA) for = 0.1, 0.05 and 0.001. When this decay is allowed, it will contribute to the 
Higgs total width thereby modifying the BRs into the other channels. In particular, it will modify the 
signal strength = r(h ^ 77)/rs'M(^ ^ 77 ): which is measured to about 10% precision (see for 
example Ref. [48]). We plot in Fig. 6. We thus see that the constraint on from the 8 TeV LHC 
is of the order of 0.01 if tra < ^/i/2. 

SVQ model: We consider a BSM extension, which we call the SVQ model, with an SU{2) singlet A, 
and one SU{2) doublet vector-like fermion 7 / = ^Pl,r = ('0iL,i?5 '02L,i?)^ with hypercharge V/,. To the 
SM Lagrangian we add 

H + '4)ilp')p - ihaA^^^iP - - ^A‘^ - ^A^H^H , (5) 

where the gauge interactions of the ip are understood and are not explicitly shown. For = 1/6 one 
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Figure 6 : BR {h AA) (left) and (right) as a function of ttia for SVU model. 


can add the terms + y'^^pLHdR + iy 2 AAqL'ipR + h.c^, which we will not consider here but later 

in Sec. 3.2.3. As in the SVU model, there are no decays to a pair of SM fermions, but unlike there, in 
this model A W~^W~ decay is also possible through the VLF loop, in addition to gg, 77 , Zj and 
ZZ modes. The expressions for the effective couplings of the A to two SM gauge-bosons are given in 
App. B. We take = 1/6 as an example. 

In Fig. 4 we plot the BR of A into 77 , Z 7 , ZZ and W^W~ modes. As in SVU model, the BR into 
gg remains almost constant at around 0.99 for ttia ^ 300 GeV. As the ' 0 i' 02 bF coupling {g) is greater 
than the 'ipi'ipiZ couplings {g/cw){T 3 — Qs^), the BR into WW is larger than into ZZ. Again, for the 
same reasons explained in the SVU model, the exclusion limits from the 8 TeV LHC in the 77 , ZZ, WW 
channels are rather weak in this model also. 

The a{gg -A- A) in this model is twice of what was obtained in the SVU model because there are two 
degenerate VLFs in the loop. The VLFs are degenerate since no Yukawa terms involving the SM Higgs 
can be written down that can split the masses after EWSB. Since no couplings to a pair of SM fermions 
exist, there are no 6 -quark initiated production processes possible. 

3.2 Models with A,H in SU(2) doublets (2HDM) 

In the 2HDM we have two scalar doublets, 4>i and 4 > 2 , which we take to have hypercharge +1/2. 
The physical neutral states are two CP-even scalars {h and H) and a CP-odd scalar (A). The Higgs 
Lagrangian is given by 

CD\DAif + \DA2f -V{^) , (6) 

where 

y($l, $ 2 ) =m\A\^l + + 2 ^ 2^2 - TO? 2 (^ 1^2 + h.c) + + A2($5$2)^ + 

A 3 ($ 1 $i)+$ 2 ) + A 4 + $ 2 )($|«'i) + y[+$ 2 )' + h.c.] . (7) 

In the limit when m ^2 = O 5 Lagrangian has a discrete Z 2 symmetry under which 4>i ^ —4>i, 
dR -A- —dR (with all other fields unchanged), if the down-type right-handed fermions couple only to the 
4>i and the up-type right-handed fermions only couple to the 4>2 so that there are no tree-level FCNCs 
(see for example Ref. [49]). Nonzero m ^2 softly breaks this Z 2 symmetry. We will not consider the hard 
Z 2 breaking terms ( 4 >| 4 >i 4 >| 4>2 + ^ 2 ^ 2 ^ 2^1 + h.c)^. There are eight free parameters in V. After we 
fix the minimum of the potential at (4>i) = (0, V\j\f2y and (4>2) = (0, n 2 /\/ 2 )^, with the constraint 
^ 1 +^ 2 =^^ = (246 GeV)^, the number of free parameters reduces to seven which we take to be m^, 
tan^d, a and mi 2 , in a notation that is common in the literature (see Ref. [50]). We 
parametrize the scalar doublets as 

~ 

with vi = vcos/3, V 2 = vsinP and tan/3 = ^ 2 /^ 1 . The physical mass eigenstates are: a heavy CP- 
even scalar H = pi cos a + p 2 Aiia, a light CP-even scalar h = —pisina + p 2 Cosa, a CP-odd scalar 

^We use the notation H = . 

^This is a natural choice since if these terms are zero to start with they will not be induced at the loop level even if the 
soft breaking terms are present. 
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A = —7]i sin P rj 2 cos f3, and charged scalars = —(f) f sin ^ (pf cos [3. All the effective couplings, 
relevant BRs and the cross sections in the 2HDM can be found in Refs. [50, 46]. The expressions of a, P 
in terms of the model parameters can be found, for example, in Ref. [19, 50]. It is these neutral scalars 
A, H that we are studying in this work. 

In some regions of parameter-space, rriA ~ i-c. their masses are within the experimental res¬ 
olution to distinguish them. If so, we must add the contributions from both A and H to any given 
channel; their sum is incoherent due to the different CP quantum-numbers. For instance, the experimen¬ 
tal invariant-mass resolution in the r+r“ channel is about 30% (see for instance Ref. [51]). Therefore, 
we consider two cases, one when tyia and tuh are within 30 % and add the contributions from the “degen¬ 
erate” A and , and another when they are split by more than 30 % and treat them separately. When 
they are degenerate, for the channel for instance, we have BR{A r+r“) ~ BR{H r+r“) 

in the so-called alignment limit (as will be defined precisely later), and we can use the constraints ob¬ 
tained in Sec. 2 if we interpret hif^gg shown there as ^^‘aqq + ^‘ngg BR^cj) tt) as BR{A 

r+r“) + BR{H r+r“). For the non-degenerate case, again one can make use of our results in Sec. 2 
to obtain constraints either for the H or A. 

We are interested in the case where the lighter CP-even scalar (h) is the observed 125 GeV Higgs 
boson. For this, the cos(;d — a) ^ 0 is the most favored region (see Fig. 18 of Ref. [17]). Only a small 
range of other values of {P — a) are allowed where the sign of the down-type coupling of the Higgs is 
reversed. For the 2HDM with exact Z 2 symmetry (i.e. 771^2 = 0), tan^d has an upper limit of 7 from 
perturbativity constraint (see Ref. [24]). We will work with a nonzero m ^2 which allows for larger values 
of tan^d (see Ref. [20]). We also assume that the “alignment limit” {13 — a = 7r/2) holds sufficiently 
accurately so that the h couplings are SM like to match with the properties of the observed 125 GeV 
state at the LHG as discussed in Ref. [16]. In this limit, the ^ WW and H ZZ decays do not give 
any significant constraints on the parameter space (see for example Ref. [47]). 

Depending on how the fermions couple to 4>i and 4>2, various types of 2HDM have been defined in 
the literature, some of which we discuss next. We start by discussing a 2HDM with only the SM fermions 
present, and follow it up with many examples of different ways of adding vector-like fermions. Although 
we do not discuss it here, it is possible in some models for the yntt coupling to be accidentally suppressed 
and consequently for the BR{(I) tt) to become small, as for example in the model discussed in Ref. [10]. 


3.2.1 Type-II 2HDM 

In the Type-H 2HDM the SM Yukawa couplings are replaced by 


C D -VdQL^idR - yuqL^ 2 UR + h.c. , 

where = icr^4>*. The Yukawa couplings of h, A to the SM fermions are given as. 


1 


C D-^ {yuhcaURUR — ydhsadRdR — yuC/diAuRUR — ydSjdiAdRdR + h.c.) 

v2 


(9) 


( 10 ) 


The i7-Yukawa couplings can be obtained from the h-Yukawa couplings by the replacements —c^ 

and Co, So,. We find the allowed regions of parameter space from the exclusion-limit on (7{gg 
(/)) X BR{(I) r+r“) presented by ATLAS [38, 51]. We focus on the channel as currently this is 

the most constraining one. We do this first in the 2HDM Type-H (2HDM-H) without the addition of 
any VLFs. 

In Fig. 7 we show the tree-level decays of A to SM fermions BR(A ^ 66 ,r+T“,tt) as a function of 
niA for various tan p for the Type-H 2HDM. The loop-level BR{A 77 , Z 7 ) in the Type-H 2HDM are 
shown in Fig. 16 by the dashed-black curves, and our results match with that of the Ref. [46]. We see 
that the BRs into 77 and Z 7 are smaller compared to that of the corresponding loop induced SM Higgs 
branching ratios even for tan P = 1 when the couplings of A to the SM fermions are equal to the Higgs 
Yukawa couplings. This is because the partial width r{h 77 , 7 Z), being dominated by the W loop, 
is larger than the partial width r(A ^ 77 , 7 Z) in which only the fermions contribute (see for example 
Fig. 2.10 of Ref. [46]). For larger tan^d the branching ratios are even smaller because of the increased 
T{A 66 ) and T{A r+r“) (recall that the Abb and Ar~^r~ couplings are proportional to tan/3). 
The discontinuity at ttia = 2mt in the BRs in Fig. 16 for tan^d = 1 is because of the onset of A ^ tt 
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Figure 7: BR{A r+r“,66) (left, middle) for tan/3 = 1,5,10,15,30 and BR{A tt) (right) for 
tan^d = 1, 5, 10, 15 in 2HDM Type-II model. The loop-level BR{A —> VV) in the Type-II 2HDM are 
shown in Fig. 16 by the dashed-black curves. 


/CAgg(2HDM-II) 



^Hgg( 2HDM—II) 



Figure 8: Contours of i^Agg (left) and i^Hgg (right) in the Type-II 2HDM. 


on-shell decay. For larger tan the discontinuity is smaller since the Ati coupling becomes smaller. The 
h ^ AA decay, possible for rriA < TUhl^^ is studied in Ref. [17] and we will not discuss it here. 

In Fig. 8, we plot contours of i^Agg and K,Hgg ia Type-II 2HDM. Using this, one can read-off the 
(j{gg (j)) at the 8 and 14 TeV LHC from Fig. 1 in Sec. 2. Using the channel constraints shown in 

Fig. 2 of Sec. 2 we obtain constraints on this model. In Fig. 9 we plot the 95% confidence level constraints 
on the m^-tan/3 plane, when only A is present (left), and for ttia = mn when both contribute (right). 
Ref. [51] has presented similar constraints in the m^-tanyd plane, but for the MSSM. 



150 200 250 300 350 400 450 500 200 300 400 500 600 

m^(GeV) rriA (GeV) 


Figure 9: For Type-II 2HDM, regions of the m^-tan/d parameter space (blue region) which is excluded 
at 95% confidence level from p ^ decay when only A is present (left) and when rriA and rriH are 

degenerate (right). 
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Figure 10: For the Type-X 2HDM, BR (74 ^ (left, right) for tan/3 = 1,5,10,15,30, and 

BR (74 ^ bb) (middle) for tan^d = 1, 5. The BR{A 77 , Z 7 ) for the Type-X 2HDM is shown in Fig. 27 
as the dashed-black curve. 


/^Agg(2HDM-X) /^Hgg(2HDM-X) 



Figure 11: For the Type-X 2HDM, contours of hZAgg (left) and i^Hgg (right). 


3.2.2 Type-X 2HDM 

In the Type-X 2HDM (2HDM-X) (see Ref. [28, 50] for a description of this model) all the SM quarks 
couple to ^2 and all the leptons couple to $ 1 . The Lagrangian for the model 2HDM-X is given by 

£ D - {ydQL^ 2 dR + yuQL^ 2 UR + yeh^ieR + h.c) + - V{^). ( 11 ) 

As a result, A coupling to the quarks and leptons are proportional to cot (3 and tan P respectively. In the 
Type-X model, since all SM quarks couple very weakly to A for large tan^d, (j{gg A) becomes very 
small for large tan^d. As a consequence there are no constraints from cr(pp A) x BR{A r+r“). 
The SM quark contribution to k,avv for 2HDM-X can obtained from that of 2HDM-II (see Ref. [46]) by 
replacing tan/3 with cot /3 in the Abb coupling. In Fig. 10 we show the tree-level BR(A ^ ^ bbRt). 

The BR{A VV) for the Type-X 2HDM is shown in Fig. 27 as the dashed-black curve. 

In Fig. 11 we plot contours of K,Agg and i^Hgg- From this one can read off (j{gg (j)) for 8 TeV and 
14 TeV LHC from Fig. 1 in Sec. 2 . The results for R^j^gg in 2 HDM-X are also applicable for Type-I 2 HDM 
as the SM quarks couple to i7, A in an identical fashion as in Type-X 2HDM. 

Next, we add various combinations of SU{2) singlet and doublet VLFs to the Type-II 2 HDM first, 
and to Type-X 2 HDM following that. Our goal is to study how VLFs affect the LHC production rate 
and decay BRs of the (j). There are eight different ways in which the 4>i and the 4>2 can couple to the 
VLFs consistent with the symmetries of the 2HDM-H namely 4>i ^ — 4>i and (Ir —dR (with all other 
fields unchanged). Among these eight models we will discuss only three representative ones that also 
capture the effects in the others. 

3.2.3 Type-II 2HDM with VLQ-SMQ Yukawa couplings 

Many models that address the hierarchy problem, such as for example composite-Higgs and little-Higgs 
models, have as an important ingredient off-diagonal couplings between a VLF and 3rd-generation SM 
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fermions. We discuss this possibility in a model-independent way by introducing, one at a time, SU(2)- 
singlet VLFs with EM charge 2/3 and —1/3. As an example, we show how the results obtained here 
apply to a little-Higgs model. 

MVU model : In what we call the MVU model, we introduce an SU(2)-singlet VLF pair (t/, denoted 
by the 4-spinor t/, with EM charge 2/3, and add to the 2HDM Type-II Lagrangian the following terms 


Cd M^tpip - (yi qL^i'ipR + h.c?j . 

After EWSB the mass terms for the EM-charge 2/3 fermions can be written as 


£ 


mass 


1 


{yuV2tLtR + yiVitL^R -h h.C.) + 


We define the mass eigenstates ^ and t 2 L,R^ foi* the EM-charge 2/3 quarks as 


tL,R = cosOY^Rtl^j^ - sin( 9 ^^t 2 L,i?, 
i^L,R = + COS 


( 12 ) 


(13) 


(14) 


The mixing angles and the mass eigenvalues can be found in App. A.l. Eor notational brevity we call 
simply as t, which we will identify with the SM top-quark. Constraints on the mixing from EWPT 
and a vector-like top decaying to Wb, Zt, Ht are studied in Refs. [5, 7, 52, 53]. Constraints from flavor 
observables are studied in Ref. [52]. 

The A couplings to the EM-charge 2/3 fermions in terms of the mass eigenstates are given by 


^ {yAtt^L^R + yAt2t2^2Lt2R + yAt2tt2LtR + yAtt2tLt2R) + h.C. 


(15) 


wherey^tt D cos/ 3 -yiC^s^ sin/3), 2/At2t2 = {VuS^s^cos l3^yiS^c^smj3),yAt^t = - {VuS l cos/3- 

yic^s^sm/3), and yAtt 2 = ~(VuC^cos/3 + //ic^c^ sin/3). The h couplings to the EM-charge 2/3 
fermions are given by. 


^ ^ {yhitthtR + yht2t2'^2Lt2R + yht2thLtR 


■ yhtt2^Lt2R) + h.C. , 


(16) 


where cos sina), ^^ 2^2 = cos sina), = {yuSLC^cosa- 

yic^s^ sin a), and yhtt 2 = {yuC^s^ cosa + ^ic^c^ sin a). We fix = 163 GeV [54] by choosing i/u ap¬ 
propriately, and show in Eig. 12 the contours of Rhtt = yhtt/yhu ^he yi-M^ plane. In the region to the 
left of the 0.99 contours, Rhtt approaches 1. The experimental constraint on r^u is 0.63 < Rhtt < 1-2 [55]. 

In Fig. 13 we show contours of i^ i^he mA-Mr^p plane for {tan/3^ i/u} = {1,1.4} and {5,1}, and 



Figure 12: For the MVU model, contours of Rhtt for tan/3 = 1 (left) and tan/3 = 5 (right) with yu 
chosen such that rrit = 163 GeV. 


also show as a function of yi for rriA = 1000 GeV, = 1250 GeV and tan/3 = 0.1 ,1, 5, 10, 15. 
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Figure 13: Contours of jyl for {tan^^, yu} = {1, 1.4} (left) and {5, 1} (middle) for the MVU 
model. as a function of yi, for uia = 1000 GeV, = 1250 GeV and tan^d = 0.1,1,5,10,15 is 

plotted on the right. 


For large tan^d, the mixing angles become small, which makes small. For Fig. 13, we fix y^ = 1.4 

so that rrit is close to its experimental value, and once a specific choice of yi is made, rrit can be fixed 
exactly by choosing yu slightly differently; the resulting change in such differences in yu is 

insignificant. 

The fermionic decay BR for rriA < + TTit) will be largely unchanged from the Type-II 2HDM 

plots shown in Fig. 7. However, if rriA > + '^t) the A t 2 t decay becomes kinematically allowed. 

In Fig. 14 we plot BR(H ^ tt), BR(H ^ 55), BR(H ^ gg) and BR(H ^ ^ 2 ^), for = 1 TeV, yi = 1 
and tan^d = {1, 5} with yu fixed such that rrit is at the physical value. BR(H ^ 77 , Z 7 ) do not change 
by much from the 2HDM-II case. 



Figure 14: Contours of BR(H ^ tt) (black), BR(H ^ 55) (blue), BR(H ^ t 2 t) (red), BR(H ^ gg) 
(green) with Mr^p = 1000 GeV, yi = 1 for tan/3 = 1 (left) and 5 (right), for the MVU model, with yu 
chosen such that rrit = 163 GeV. 

As an example, we apply these results to a concrete model that stabilizes the electroweak scale, has 
a 2HDM structure, and has vector-like fermions, namely the SU(6)/Sp(6) little-Higgs model by Low, 
Skiba and Smith (LSS, Ref. [11]), which we analyze in detail in Ref. [10]. Among the various sample 
points that are listed in App. B in Ref. [10] that satisfy all constraints including precision electroweak, 
we consider here the sample-points 1 and 2. For the sample-point-1, the two lightest VLFs are the t 2 
with a mass of 1218 GeV, the 52 with a mass of 1315 GeV, and we have tan/3 = 1.36, ijia = 1671 GeV, 
yi = 1.7, yu = 1.2 and rrit ~ 164 GeV."^ Keeping only the lighter t 2 since the ts is somewhat heavier, a 
good approximation is obtained by considering the addition of only a singlet EM charge +2/3 state 7 / 
as introduced in Eq. (12). Ignoring the smaller 52 contribution, the /^X^^ due to the ^2 can be read off 
from the tan^d = 1 curve of the rightmost panel of Eig. 13 to be approximately 0.4. This is about 10% 
of the SM-fermion contribution. 

‘^From Ref. [ 10 ] we have the following for the LSS model: there we had tan /3 = vxjv^ while in this paper we have 
tan (3 = V2IV1, and therefore tan (3 here is related to that of Ref. [ 10 ] via tan (3 = (1/tan ( 3 ^^^); yi is given by yi = yf‘^‘^023, 
and for point-1, since yi ^ y4, to a very good approximation S14 1 and C14 0; also mt ~ C 23 y 2 V 2 lV^ in the limit 

where ts is decoupled away, i.e. yu = y2C23, and C23 ~ 0 . 9 . The 62 is an SU( 2 )-singlet since it does not mix with the other 
charge - 1/3 states. 
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MVP model : In the MVD model, we introduce an SU(2)-singlet VLF pair denoted by the 

4-spinor x, with EM charge -1/3, and add to the 2HDM Type-II Lagrangian the following terms 

Ca = M^xX- {y 2 Ql^iXr + h.c.) . (17) 

The mass eigenstates, 6 ^ ^ and h 2 L,R for the EM-charge —1/3 fermions are defined in the same way as 
in Eq. (14) with the mixing angles, 0 ^The mixing angles and the mass eigenvalues can be found in 
App. A.l. The A couplings to the EM-charge —1/3 fermions are obtained in a similar way as in Eq. (15), 
with the replacements ^c^sin/3, yi y 2 - Similarly, the h couplings to the EM-charge —1/3 

fermions are obtained from Eq. (16), with the replacements yu cos a ^ —yd sin a and yi ^ ^ 2 - As in the 
case of charge 2/3 fermions, we choose yd such that = 4.2 GeV [56]; yhbb stays close it’s SM value. 

In Eig. 15 we plot contours of ly\ in the tua — plane for {tan^d, yd} = {1, 0.03} and {5, 0.12}, 
and as a function of y 2 for uia = 1500 GeV, = 1000 GeV for tan^d = 1 , 5, 10 , 15. 




m^(GeV) m^(GeV) y2 


Eigure 15: Gontours of /y^ for {tan/3, yd} = {1, 0.03} (left) and {5, 0 . 12 } (middle) for the MVD 
model. with y 2 for ttia = 1500 GeV, = 1000 GeV and tan/3 = 1 , 5, 10 , 15 is plotted on the 

right. 

As an example we consider again the LSS model, but now the sample-point -2 in App. B of Ref. [10], 
with the lightest VLE being the 62 with a mass of 947.5 GeV. The 62 is an SU(2)-singlet state and 
does not mix with the other charge -1/3 states. Eor this point, ttia = 1671 GeV, tan/3 = 1.36, y 2 = 
1.422, C 23 = 1.15. RAgg cau be read off from the tan^d = 1 curve of the rightmost panel of Eig. 15 to be 
approximately 0.3. 

MVQ model: Eor the MVQ model, we add an SU( 2 ) doublet VLE pair denoted by the 

4-spinor Q', and add to the Type-II 2HDM Lagrangian the terms 

C D Mqq Q'Q' + (^MqQ QlQ'r - yiQL^2tR - y2Q'L^ibR + h.c.^ . (18) 

In the following we show only the top-sector since this is usually the dominant piece in BSM models, and 
we therefore suppress the bottom-sector. At the outset, we diagonalize the VLE masses by redefining 
the Q and Q' fields by an orthogonal rotation to get an equivalent Lagrangian given by 

C D Q'Q' + {-yfqMR - ri^Q'L^ 2 tR + h.c.) , (19) 

where we show the 2HDM top Yukawa coupling also since its effective coupling is now changed, with 
^ ^(M2 q + M2q), yf = {y^MQQ-y,M,Q)/M^Q% yf ^ {y^M,Q + y,MQQ)/ which imply 

UtF = iVu - yiMqo/M qq )/+ (MqQ/MQQY and yl^ = y/f M^q/Mqq + yi 

The Rcj^gg due to the t'^h' in the MVQ model are qualitatively similar to the MVU case presented 
earlier. As an example, let us consider again the LSS model sample-point -1 in App. B of Ref. [10], for 
which we have = 0, y^^ ^1.3 and y}^ ~ 0.5, which gives RAgg ~ 0.03. The doublet-VLQ contribution 
in this case is thus very small compared to the SMQ contribution. 
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3.2.4 Type-II 2HDM with VLQ-VLQ Yukawa couplings 

Here, we add SU(2) doublet and singlet VLFs with SM-like hypercharge assignments, and write Yukawa 
couplings between them both involving the 4>i^2- Although there could be Yukawa couplings between a 
VLF and an SMF also present, we do not write them here for simplicity, and their effects are investigated 
separately in Sec. 3.2.3. 

MVQDll niodel: To the Type-II 2HDM we introduce one doublet VLQ, 'ip = (V^i, '^ 2 ) with hypercharge 
and one singlet VLQ (y) with hypercharge (V/; ~ 1/2) so that VLF couplings with 4>i are allowed. 
The additional Lagrangian terms to the 2HDM-II are 

C D xi)ilpip + xilpx - (yi^L^iXi? + Vii’R^iXL + h.c) - - M^XX- (20) 


We can also write the terms 'ipL^2XR ^^^<1 '^r^ 2 Xl^ which we do not add here but will consider them 
subsequently as another model. These terms are forbidden if y ^ —y under the Z 2 symmetry of 
2HDM-IL The terms involving h, A and VLFs after EWSB are 


- 1 - - V - - 

Cd- - M^XX + -^Asmp{iyitp2LXR + iyii’2RXL + h.c.) - — cos p{yi'ip2LXR + yi^2RXL + h.c.) 


+ -^hsma{yii;2LXR + yii’2RXL + h.c.) - -Xh cos a{yx'ip2LXR + yii’2RXL + h.c.) . 


( 21 ) 


Gauge interactions of the VLFs are present and not shown explicitly. 7/^2 and y mix after EWSB, while 
pji is itself a mass eigenstate. We define the mass eigenstates Q and (2 as 

^2L,R = ClL,R COS 0l,R “ C2L,R An 0l,r, (22) 

Xl,r = CiL,i? sin + C2L,R cos Ol^r, (23) 


where the mixing angles Or and Or are defined in App. A.2. In terms of these mass eigenstates, the 
Lagrangian in Eq. (21) can be written as 

£ D - yfjiiACiRjR + h.c) - MiOO - M^Ripi + + eQiAf^QxnCi 

- yij {hCiiCjR + h.c) - y^j {hCiiCjR + h.c) , ( 24 ) 

where i,j = 1,2 and yf^'s are given in App. A.2. We take the yi and yi to be real, enforcing CP invariance 
in the BSM sector. The relative sign between yi and yi in Eq. (20) is physical for the following reason. 
If we want to get rid of this relative sign we need to make the transformations xl ~Xl and xr XRi 
or Xl ^ Xl and xr ~Xr- la either case, the changes its sign and is therefore a physical effect. 
Eor chiral fermions, the sign of the mass term is not physical since one can rotate it away by the above 
transformations. 

Instead of the y (with hypercharge (V/; —1/2)), if we consider a VLE (say ^) of hyper charge (V/; + l/2), 
we get a different model where the ^ couples to the instead of the 4>i. This model will have similar 
phenomenology as MVQDn model, which we discuss later. 

The effective couplings for this model are given in App. B. When yi = ^i, in addition to CP 
invariance, the Lagrangian in Eq. (21) is also invariant under P and C individually, with A transforming 

PC ~ 

as A —A, A —)► —A. This implies that the VLE contribution to Raw is zero since AV^^V^^ is not P 

invariant (although it is CP invariant). Also, the VLE contributions are maximum for when 

the mixing between the VLEs {'ip 2 and y) is maximum. We will take Mr^p and to be equal from now 


In Eig. 16, we plot BR(A ^ W), for V/; = 1/6 as an example, which is the SM quark-doublet 
hypercharge assignment. The tree-level decays to SM fermions BR(A -> 66 ,r+r“,tt) are unchanged 
from what is shown in Eig. 7 for the Type-II 2HDM. We see that for small values of tan/3 the VLE 
contribution to BR(A ^ VV) is small compared to the 2 HDM-IL This is because yij^s are proportional 
to sin /3. Eor large tan ^ and for large m^, the VLE contributions to the BR(A ^ 77 ) become significant. 

In Eig. 17, we plot contours of raqq for = 800 GeV, 1700 GeV. Eor comparison we have also 
plotted the corresponding contours in 2HDM-IL Using this, one can read-off the a{gg A) at the 
8 and 14 TeV LHG from Eig. 1 in Sec. 2 . Eor comparison, the corresponding contours in the Type-II 
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Figure 16: BR (74 ^ 77 ) (top-panel) and BR (74 ^ Z 7 ) (bottom-panel) with = 1000 (GeV), 

tan/3 = 1 (left) and 30 (right) in MVQDn model (solid-black) and in 2HDM Type-II (dashed-black). 
BR{A //) in Type-II 2HDM are as shown in Fig. 7. 
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Figure 17: Contours of K.Agg for = 800 GeV (left), 1700 GeV (middle), yi = 0.5, = 1 for 

MVQDii model. 
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Figure 18: For the MVQDn model, yfi (red), i/ii (blue) as a function of tan^d (left); K,Agg as a function 
of tan/3 for ttia = 300 GeV (middle) and 600 GeV (right), with yi = 0.5, yi = 1 and = 800 GeV 
(blue), 1000 GeV (green) . 
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Figure 19: Contours of i^Hgg for yi = 0.5, yi = 1, for = 800 GeV (left), 1700 GeV (right) for 

MVQDii model. The corresponding contours in Type-II 2HDM are shown in Fig. 8. 


2HDM (without the VLFs) are shown in Fig. 8. In Fig. 18 (left) we plot y^^ and yf^ (defined in Eq. (24)) 
in the alignment limit (^d — a = 7r/2), which shows that the h couplings to the VLFs become very small 
as tan/3 increases. Thus, the VLFs can modify (j{gg A) and T{A VV) significantly, while the h 
remains SM-like as required by the LHC measurements of the 125 GeV state. We find that the VLF 
contributions partially cancel the SM fermion contributions for a range of low tan /3 values and for some 
ranges of the m^, while for larger tan P the effective couplings always increase compared to the 2HDM-IL 
To illustrate this point more explicitly, we plot i^Agg as a function of tan^d in Fig. 18 for tua = 300 GeV 
and 600 GeV. The constraint on the 2HDM was nontrivial only for large tan P (see Fig. 9). Therefore, for 
large tan^d, since the i^Agg is bigger for this model compared to 2HDM (see Fig. 18), and the tree-level 
BR from which the tightest constraint appears is almost unchanged, the constraint on this model 
will be tighter. In Fig. 19, we plot contours of i^Hgg foi* , in the alignment limit. Corresponding 

contours in Type-II 2HDM are shown in Fig. 8. From this, one can also obtain a{gg -A- H) from Fig. 1. 

MVQU 22 model: We introduce one doublet VLQ (-0) with hypercharge V/, and one singlet VLQ (^) 
with hypercharge V/; + 1/2, which couples only to 4>2. We add the following terms to the 2HDM-II 
Lagrangian 

C D ipupip + - y2'fpL^2^R - y2'4’R^2^L + h.c - (25) 

Here we do not include the terms and as their effects have been considered in MVQDn 

model. As the BR(A ^ W)s do not change much compared to the 2HDM-H case, we do not show 
them here. Instead of the ^ (with hypercharge V/; + 1/2) if we consider a VLF (say x) of hypercharge 
— 1/2) we get a different model where the y couples to the 4>2 instead of the 4>2. This will give 
similar effects to what we consider here. 

Similar to MVQDn model we diagonalize the mass matrix by an orthogonal rotation and define the 
couplings yfj. The mass eigenvalues, mixing angles and yfj^s for this model can be found in App. A.2. 
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^Agg(yi=0-5, = A4=800) ^Agg(ji=0-5, ji = l,M^=1700) 



Figure 20: Contours of i^Agg for ^2 = 0.5, ^2 = 1, for 
MVQU 22 model. 


^Hgg(yi=0.5, M^=800) 



M^ = M^ = 800 GeV (left), 1700 GeV (right) for 
^Hgg(3'^i~^‘5) Jh~F A^^=1700) 



Figure 21: Contours of ^Hgg for 1/2 = 0.5, ^2 = 1, for = 800 GeV (left), 1700 GeV (right) for 

MVQU 22 model. 


The effective couplings for this model are given in App. B. As in MVQDn model, the k,avv becomes 
zero when ^2 = ^ 2 - In Fig- 20 we plot contours of K.Agg m m^-tan/^ plane. In the MVQU 22 model 
the VLF contributions to i^Agg are very small for yi = 0.5 and = 1, and therefore we do not show 
it explicitly. This is particularly so for large tan/3 because ^^^’s are proportional to cos^d which become 
small as tan/3 increases. Similar conclusions hold for i^Hgg- In Fig. 21 we plot i^Hgg using which one 
can read-off the cr^gg H) from Fig. 1 by reading K.Agg there as t^Hgg as mentioned earlier. Since i^Agg 
and i^Hgg do not change much compared to the 2HDM-II, constraints on the m^-tan /3 plane will almost 
remain same as in the 2HDM-II case. Thus, VLFs if realized as in MVQU 22 model have little impact 
on the observables we consider here. 

MVQU12 model: We introduce one doublet VLQ (-0) with hypercharge and one singlet VLQ (^) 
with hypercharge (V/; + 1/2)- We consider the case where couples only to <I>i and couples only to 
^ 2 - To the 2HDM-II Lagrangian, we add 

C D + h.c) - (26) 

We get different models if instead of the couplings above, the ipR couples to 4>i and couples to ^> 2 , 
or, if instead of ^ we introduce a VLF singlet (say y) with hypercharge {Y^ ~ 1/2) with couplings to <I>i 
and ^ 2 - All these models have similar phenomenology as MVQU 12 model. 

The mass eigenvalues, mixing angles and yfj^s for this model can be found in App. A.3. The effective 
couplings for this model are given in App. B. In this model, the effective couplings do not reduce to zero 
for yi = ^ 1 , unlike in MVQDn and MVQU 22 models, as there are no additional P and C symmetries 
in the VLF sector. In Fig. 22, we plot the BR(A ^ VV), BR(A ^ 66, , it) for an example choice of 
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Figure 22 : BR (74 ^ 77 ) (left) and BR(A ^ Zj) (right) with = 1000 GeV for tan^d = 30 

for MVQU 12 model (solid-black), and the corresponding variation in the Type-II 2 HDM (dashed-black). 
These BR for tan (3 = 1 and the BR{A rr, 66 , tt) are not explicitly shown here as they are identical 
to those in Figs. 16 and 7 respectively. 

^Agg(ji = B Jy = l,M^=800) KA^(yi = l, y^ = l, M^=1700) 




Figure 23: Contours of K.Agg foi* Hi = 1, ^1 = 1, foi* ^7 = ~ (left) and 1700 GeV (right) 

for MVQU 12 model. 


Yr^^J = 1/6. The BR{A 77 , Zj) for tan = 1, = 0.5, = 1 and the tree-level BR{A rr, 66 , tt) are 

not explicitly shown in Fig. 22 as they are identical to those shown for the MVQDn model in Fig. 16, 
and Type-II 2HDM in Fig. 7 respectively. In Fig. 23 we plot contours of i^Agg for ^1 = ^1 = 1 and 
= 800 GeV and 1700 GeV. From this, one can obtain (j{gg A) at the 8 and 14 TeV LHC 
from Fig. 1 in Sec. 2. For low values of tan/3 the effective coupling increases compared to the 2HDM-II 
case, while for larger values of tan /3 the effective coupling decreases compared to the 2HDM-II. To show 
this more explicitly, we plot i^Agg with tan^d in Fig. 24. The decreased coupling is due to a destructive 
interference between the contributions from SM fermions and the VLFs. If we reverse the sign of yi or 
^ 1 , we get the opposite effect; for low values of tan/3 the effective coupling decreases compared to the 
2HDM-II while for larger values of tan/3 the effective coupling increases compared to the 2HDM-II. In 
Fig. 25 we plot contours of tZHgg in the alignment limit. From this, one can also obtain a{gg H) from 
Fig. 1 by reading ^Agg there as tZHgg as mentioned earlier. 

In Fig. 26 we plot the region of the my^-tan^d parameter-space which is excluded at 95% confidence 
level for two cases, when only A is present, and when A and H are degenerate and both present. 
For comparison, we have also plotted the corresponding limit for the 2HDM-II case. We see that the 
constraints are loosened compared to the 2HDM-II due to the presence of VLFs. This happens because 
of the reduction of hZAgg {i^Hgg) compared to the 2HDM-IL 

Next, we add VLFs to the Type-X 2HDM and study the phenomenology of the neutral scalars. 

3.2.5 Type-X 2HDM with VLQ-VLQ Yukawa couplings 

MVQDXii model: To the 2HDM Type-X model in Eq. (11), we introduce VLFs in a similar fashion as 
in MVQDii model as a representative case, and call it MVQDXn model. The other ways of coupling 
VLFs similar to MVQU 22 or MVQU 12 model will be qualitatively similar to our results here. We 
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Figure 24: tvAgg with tan^^ for ttia = 300 GeV (left) and 600 GeV (right) with = 1, = 1 and 

= 800 GeV (blue), 1000 GeV (green) for MVQU 12 model and 2HDM-II (dashed-black). 




^Hgg(34~l> 800) 



^Hgg(34 = F yi = ^> M^=1700) 



Figure 25: Contours of i^Hgg for Hi = 1, ^1 = 1, for = 800 GeV (left), 1700 GeV (right) for 

MVQU 12 model. 
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Figure 26: For MVQU 12 model, regions of the m^-tan^d parameter-space excluded at the 95% CL from 
r+r“ decay when only A is present (left), and when A and H are degenerate and both present 
(right), with yi = yi = 1, = 800 GeV (dark-blue region), 1000 GeV (light-blue and dark-blue 

regions). All shaded regions are excluded in the 2HDM-II. 
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tanyS=30, 3 / 1 = 0 .5, 3 / 1=1 
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Figure 27: BR (74 ^ 77 ,^ 7 ) with = 1000 GeV (solid-black) for tan^d = 30 for the 

MVQDXii model, and the corresponding variation in the 2HDM-X (dashed-black). 
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Figure 28: Contours of tvAgg for yi = 0.5, yi = 1 , for = 800 GeV (left), 1700 GeV (right) for 

MVQDXii model. The corresponding contours in Type-X 2HDM is shown in Fig. 11. 


introduce a doublet VLQ ip = (' 0 i,' 02 ) with hypercharge V/;, and a singlet VLQ (y) with hypercharge 
(V/; — 1/2) which couples only to ^ 1 . To the 2HDM-X Lagrangian we add 

£ D Ipilp^ + X0X - {vii’L^iXR + yii’R^iXL + h.c) - - Max- (27) 

The effective couplings of A with VLFs are same as in MVQDn model and can be read off from App. B. 
In Fig. 27 we show BR(A ^ VV) including the VLF contributions for the MVQDXn model, and the 
tree-level BR(A ^ r+r“, 66 ,tt) is unchanged from what are shown in Fig. 10. BR(A 77 ,^ 7 ) for 
tan^d = 1, = 0.5, yi = I are almost identical to the 2HDM values shown in Fig. 16 and are therefore 

not shown explicitly in Fig. 27. For tan^d = 30, BR(A ^ 77 ,^ 7 ) is increased compared to 2HDM-II, 
because for large tan/3, r(A ^ bb) becomes much smaller in 2HDM-X. 

In Fig. 28 we plot contours of i^Agg- The hif^gg contours in 2HDM-X (without VLFs) are shown in 
Fig. 11 . Using these plots, one can read off cr{gg A) for 8 TeV and 14 TeV LHC from Fig. 1 in Sec. 2 . 
As expected, for large tan/3, tZAgg is significantly larger in this model compared to 2 HDM-X since the 
VLFs contribute substantially while the SM quark contributions alone are very small. In order to show 
explicitly how large the change is, we plot K,Agg as a function of tan (3 for rriA = 300 GeV and 600 GeV 
in Fig. 29. In Fig. 30 we plot contours of hZHgg in m^-tanyd plane in the alignment limit. From this, one 
can also obtain a{gg H) from Fig. 1. 

3.2.6 Type-II 2HDM with VLL-VLL Yukawa couplings 

MVLEii model: Vector-like leptons do not contribute in gg A, but can contribute in A ^ 77 , Z 7 . 
We show the effect of VLLs in a simple model similar to MVQDn model, but with VLLs instead of VLQs. 
We introduce one doublet VLL ( 7 /^) with hypercharge, V/; and one singlet VLL (y) with hypercharge, 
V/; — 1/2. The Lagrangian we consider is exactly the same as in Eq. (20), except here the VLLs ^ and 
y do not couple to gluons. The effective couplings are the same as for MVQDn model except for color 
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Figure 29: i^Agg with tan/3 for rriA = 300 GeV (left) and 600 GeV (right) with yi = 0.5, yi = 1 and 
= 800 GeV (blue), 1000 GeV (green) for MVQDXn model and 2 HDM-X (dashed-black). 

^Hgg(yi=0.5, J/i = l,M^=800) /^Hgg(yi=0.5, 3>i = 1,M^=1700) 




Figure 30: Gontours of K,Hgg for yi = 0.5, yi = 1 , for = 800 GeV (left), 1700 GeV (right) for 

MVQDXii model. The corresponding contours in Type-X 2HDM is shown in Fig. 11. 


factors. As an example, we choose V/; = —1/2 and plot BR(A ^ 77 ) as a function of rriA in Fig. 31, 
with = M^= 500 GeV, for tan/3 = 1 and 30. We see that the effect of VLLs is qualitatively similar 
to vector-like quarks; for low tan P the effect of VLLs is negligible while for large tan P and large rriA 
VLL contributions are significant. Near ttia = 1000 GeV, the VLL contribution is quite large due to 
them going onshell for our choice of VLL mass of 500 GeV. BR(A ^ Z 7 ) will show the same behavior. 

4 Conclusions 

Many theories beyond the standard model (BSM) contain new CP-odd and CP-even neutral scalars 
(j) = {A^H} and new vector-like fermions {'ipyh)- We study the LHC phenomenology of p taking into 
account 2pvL contributions to pgg^ ^77 and couplings at the one-loop level. 


tany0=l, yi=0.5, yi=l 


tan/?=30, yi=0.5, yi=l 




Figure 31: BR(A ^ 77 ) with = 500 GeV, yi = 0.5, /7i = 1 for tan^d = 1 (left) and tan/3 = 30 

(right) for MVLEn. 
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In Sec. 2 we write an effective Lagrangian with (j) coupled to standard model (SM) gauge-bosons and 
fermions. We focus only on (p Yukawa couplings to third generation SM fermions, namely since 

these are usually the bigger couplings in most BSM extensions. The couplings of the A to standard model 
W, Z gauge bosons (i.e. AVV couplings) cannot occur from renormalizable operators in a CP-conserving 
sector, but can be induced as loop-generated non-renormalizable operators. These operators are induced 
by SM fermions and also the heavy 'ipvL- In Sec. 2 we present model-independent results that are useful 
whatever be the origin of these effective couplings. In Fig. 1 we present the 8 TeV and 14 TeV LHC 
gg ^ (p (gluon-fusion channel) cross-sections as a function of the effective couplings. We also obtain 
limits on the effective couplings from the 8 TeV LHC data on the 77 , r+r“ and tt modes. We do not 
include the bb decay mode and the 6 -quark associated production channels in this work. 

We define some simple models in Sec. 3 that are representative of BSM constructions as far as 
the phenomenology of p is concerned. These models include p and pvL in the singlet and doublet 
representations of SU{2). In the doublet case, we focus on the two-Higgs-doublet (2HDM) Type-H and 
Type-X models. We compute the pgg^ ^77 and pZ^ effective couplings induced by the SM fermions and 
vector-like fermions at the 1-loop level and present analytical expressions for them in App. A. For the 
various models we define, we present the effective couplings ^^ 077 , and, BR{A 77 ,^ 7 ) 

and BR{A //) for / = {r, 6 , t} as a function of the model parameters. From the hif^gg and the BR 
into one of these modes, one can see if a point in parameter-space in a given model is allowed by the 
8 TeV data from our plots in Sec. 2 . One can also read-off the gluon-fusion cross-section at the 8 TeV 
and 14 TeV LHC from Fig 1 . Interestingly, for some of the 2 HDM cases we studied, we find that the 
addition of vector-like fermions loosens the constraint compared to the 2HDM alone, and allows more 
of the parameter-space. This can be seen for instance in Fig. 26. The 14 TeV LHC gluon-fusion c.s. of 
the p and its BRs in the different models we present should be useful in identifying allowed regions of 
parameter space and promising discovery channels of the p. In this context, it is interesting to explore 
the possibility of the p being the state responsible for the recent 750 GeV excess of diphoton events at 
the 13 TeV LHC [57, 58]; this is the subject of a separate work. 

Acknowledgments: We thank Gautam Bhattacharya and Dipankar Das for a discussion on the unitarity 
bounds in 2 HDM, and V. Ravindran for a discussion on bb fusion. 


A Couplings, masses and mixing angles in various models 

Here we provide the explicit expressions for the mixing angles and the mass eigenvalues in the different 
models we have defined in Sec. 3. We also provide expressions for the /^ij’s, yij^s defined in Eq. (24). 
Sec. A.l contains the mixing angles and the mass eigenvalues for MVU^ MVD and MVQ models. 
Sec. A .2 contains explicit expressions for ^^j’s and for MVQDn, MVQU 22 and MVQDXn models. 
Sec. A.3 contains explicit expressions for i/ij^s and Kij^s for MVQU 12 model. In what follows we will use 
the notations cl^r = cos Or,sl,r = sinC /3 = cos^d and = sin/3. 


A.l MVU, MVD, MVQ models 

In this section we give the mixing angles and the mass eigenvalues for MVU^ MVD and MVQ models. 
The mixing angles for MVU model are given by. 


tan26>£^ = 


2V2yiViM^ 


r, tan26>^ = 


2V2yiyuViV2 


vlvl - 2M2 + yRl ’ “““ " ^ vlvl - 2M2 - yRl' 
The mass eigenvalues for the EM charge-2/3 fermions in MVU model are given by. 


= 2 





yu 


^V2 + M^] ^^-±viT^^^V2-M^] 




(A.28) 


(A.29) 


The mixing angle and mass eigenvalues for MVD model are obtained from Eq. (A.28) and Eq. (A.29) 
by the replacements yi y 2 ^ ^u '^2 ^ and M^ M^. The mixing angles 0^ for MVQ model 
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are given by, 


tan 2^^ = 


2V2yf^V2M^^^ 


tan 2^^ = 
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The mass eigenvalues for the EM charge-2/3 fermions in MVQ model are given by. 
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A.2 MVQDn, MVQU 22 models 

In this section we give the expressions for the yij^s and Kij^s for MVQDn, MVQU 22 models. The 
couplings K,ij defined in Eq. (24) for MVQDn, MVQU 22 models and also for MVQDXn model 
are given by, ^n = (^/ch^)[(T^/2)(c| + c\) - Qs‘^], K 22 = {g/cw)[{T^/‘ 2 ){sl + - Qs‘^], 1^12 = 

— {g/cw){T^/‘^){sLCL + srCr). The mass eigenvalues Mi ^2 (in Eq. 24) for MVQDn model are given by. 


^ 1,2 — 


\ \J{M4, + M ^)2 + ^c^v'^iyi - yif ± 


^ (M^ - M^f + ^v‘^cl{yi + 2/1)2 


and the mixing angles Or^r for MVQDn model are given by. 


(A.30) 


tan26>L = 


2\/2vcis{yiM^ + yiM^) 

2{Ml - M2) - v2s2 (^2 _ yT^ ’ 


tan 20 R = 


2\/2vci3{yiM^ +yiM^) 

2(M2 - M2) + v^sl{yl - yl)' 


(A.31) 


The mass eigenvalues and the mixing angles for MVQU 22 model can be obtained from Eq. (A.30) 
and Eq. (A.31) by the replacements yi y 2 and s^. The couplings ^^’s (in Eq. 24) for 

MVQDn model are given by, y^-^ = { 1 /\/ 2 )s^{-yiCLSR + yiSRCR), = { 1 /V^)si 3 {yiSLCR - yiCLSR), 
yi 2 = -{'^/V^)si 3 {yiCLCR^yiSLSR), yn = Q/V^)sf 3 {yiSLSR^yiCLCR). T\ieyfj's\nMVQU 22 model can 
be obtained from the ^^’s in MVQDn model by the replacements yi y 2 and cp. The couplings 
Vij (ill Eq. 24) are given by, y^^ = -{1/V^)sa{yiCLSR ^yiSRCR), y ^2 = {^/V^)sa{yiSLCR ^ yiCRSR), 
yi 2 = -{^/V^)sa{yiCLCR-yiSLCR), y^i = -{l/V^)sa{-yiSLSR^yiCLCR). The^f^’sinME(5[/22 model 
can be obtained from ^^’s in MVQDn model by the replacements yi y 2 and So^ The cou¬ 

plings y^j (in Eq. 24) can be obtained from the ^^ ’s in MVQDn model by the replacements, So^ —c^ 
in case of MVQDn and 5 q, ^ —5q, for MVQU 22 model. 


A.3 MVQU 12 model 

In this section we give the expressions for the yij^s and for the MVQU 12 model. The couplings nij 
for MVQU 12 model are same as in MVQDn model. The mass eigenvalues are given by 


Mi^2 = + -v‘^{yici3 - ± J (M^ - M^)^ + -v‘^{yici3 + yiS/s) 


(A.32) 


and the mixing angles Or^r are given by 


tan26>L 


2y/2v{yicyM^ + yispM^) 

2(M2-M2)-t,2(y2g2 _y2^2)> 


tan 20 R 


2V2v{yiC/3M^ + yiSpM^) 

2(M2 - m 2) + t> 2 (y 252 _ y2g2) ' 


(A.33) 


The couplings yfj are given by, = {l/V^){yiSi 3 CLSR + yiCjsSLCR), = -(l/\/2)(yiS/}SLCi? + 
yiCpCiSR), 2/^2 = {^/V^){yiSpCLCR — yiCySRSn), y^i = —{^/'/^){yiS 0 SLSR — yiCpCLCR). The cou¬ 
plings y'l^ are given by, y5*i = (l/\/2)(-yis„CLSK -h yiCaSLCR), y ^2 = (l/\/2)(yiSaSLCi{ - yiCaCLSR), 
Vu = -{^/V2){yiSaCLCR + yiCaSLSR), = {l/V^){yiSaSLSR + yiCaCLCR). The y^n can be obtained 
from y^ ’s by the replacements s„ ^ —Ca and Ca sa- 
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B The K,(i)gg^(j)j^^fj)Z'y,Aww,Azz effective couplings in various models 


In this section we give the expressions for the K,{(f)gg^ ^ 77 , 0 ^ 7 , Azz, Aww} Ihe various models we have 
considered in Sec. 3. 


i^cfygg^ The 1-loop cxprcssious for the (j)gg and ^77 amplitudes and respectively, with 

4> = {h, A} are given here. These amplitudes are induced by quarks whose effective Lagrangian can 
be written as D f f ^y(f)f f4^f f ‘ Defining Vf = mj/m^ and with / running over all colored fermion 
species with mass rrif and Yukawa couplings ^ 0 //, and with the electric charge of the fermion (/) denoted 
by Q/, the general expressions for hCf^gg and are given as 

= 2e (B.34) 

f ruf f 

with F.';>{r,)=4r;(^ dy’ 

and g{x^y) = (1 — Axy) for the CP-even scalars (h^H) and 1 for the CP-odd scalar (A). Here M is a 
mass scale defined in Eq. ( 1 ), which we set to 1 TeV for numerical results. Compared to tiff^gg has 

an extra factor of 1/2 which compensates for our definition of r(0 ^ gg) in Eq. (2) with a relative factor 
of eight compared to r(0 ^ 77 ) while the actual color factor is really two. The expressions for in 
Eq. (B.34) match with the closed form expressions given in Ref. [59]. 

KdAZj' Here we give the general expressions for hZAZj (defined in Eq. 1 ) for the different models we have 
considered. Eor SVU and SVQ models, 
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^ v(2)/ 


k.aZ'i = 2e—^NlQi{Tl - Qis\r)yA — Fl%{ri,rz) , 

C\Y . TTli 


nl nl — y 

with F[^l{ri,rz) = ^Vi J ^y J 


1 


n + {rz - l)xy + rz{x‘^ - x) 


Eor SVU model only one VLE contributes to /^az 7 - Foi* MVQDh, MVQU 22 , MVQU 12 and MVQDXn 
models kazj = i^azj + i^azj + ‘^azj + 


^AZ 


1 


n 


ij 

AZ-i 


2cN qQ ii^aya 


M 

rrii 






x) + TjX -h {rz — F)xy + rz{x‘^ 



The couplings yij for each of the four cases are given in App. A. 2 , A.3. The expression for d /1 is a 
generalization to vector-like fermions of the expression given in Ref. [59]. 

f^Azz^ f^AWW- Here we provide the expressions for kazz and kaww for SVU and SVQ models. Eor 
the SVU and SVQ models. 


Eor 


= 2 ( A ) Ed(d - QiS^W?yA^F[%{Ti,Tz), 

\cw J . rxii 

I 

where = [(;7)e ^ 

SVQ model, kaww = 2 (^) For SVU model kaww is zero. 
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